Self-consistent field theory of polarized BEC: dispersion of collective excitation 
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We propose the construction of a set of quantum hydrodynamics equations for the Bose-Einstein con- 
densate (BEC) where atoms have electric dipole moment. The contribution of the dipole-dipole interac- 
tions (DDI) to the Euler equation is estimated. Quantum equations for the evolution of medium polariza- 
tion are constructed. The mathematical method we developed allows studying the effects of interactions 
on the evolution of polarization. The developed method may be applied to various physical systems in 
which dynamics is affected by DDI. Derivation of Gross-Pitaevskii equation for polarized particles is dis- 
cussed using point of view of quantum hydrodynamics. We shown that Gross-Pitaevskii equation appear 
at condition there all dipoles has the same direction which does not change in time. Comparison of equa- 
tion of electric dipole interaction with equation of magnetization evolution is described. A problem of 
elementary excitations in BEC, either affected or not affected by the uniform external electric field, is 
considered using our method. We show that the evolution of polarization in BEC leads to the formation 
of a novel type of elementary excitations. Detailed description of dispersion of elementary excitations is 
presented. We consider also the process of wave generation in polarized BEC by means of monoenergetic 
beam of neutral polarized particles. We compute the possibilities of the generation of Bogoliubov's modes 
and polarization modes by the dipole beam. 



I. INTRODUCTION 

After obtaining the Bose-Einstein condensate (BEC) in 
experiment with pairs of alkaline metal atoms, theoretical 
and experimental investigation of linear waves and nonlin- 
ear structures in it have been performed. In recent years 
the interest to polarized BEC has been increasing. It con- 
nected with recent experimental progress in cooling of po- 
larized atoms and molecules. At present day the BEC with 
magnetic polarization is realized on atoms 52 Cr. There are 
a lot of attempt of experimental obtaining of electrically 
polarized BEC (see review 10]). For this aim BEC 
of molecules has been cooled. Particular interest to electri- 
cally polarized BEC brought because of large dipole-dipole 
scattering length, and thereby because of both the big mag- 
nitude and large distance of interaction compare with anal- 
ogous quantities of the magnetized BEC. 

Many processes in quantum systems are determined by 
the dynamics and dispersion of elementary excitations (EE) 
0]. A law of dispersion of EE in degenerate dilute Bose 
gas was obtained by Bogoliubov in 1947 lalT-D- Many au- 
thors studied the change of the Bogoliubov's mode which 
arises when SRI interactions accounting more carefully 
- fMll . or geometry of the system is complex lfl5ll . lflr3l . 
J. In articles lfl8ll - ll22ll authors studied the influence 
of electric dipole moment (EDM) dynamics on dispersion 
of EE in BEC. The contribution of polarization in disper- 
sion law of Bogoliubov's mode was obtained in fTill - Ir22ll . 
Instability of the Bogoliubov's mode in 3D dipoles BEC 
(DBEC) with repulsive short range interaction (SRI) were 
shown in Ref. OH- HI. U. R. Fisher HI obtained 
that Bogoliubov's mode in 2D DBEC is stable for a wide 
range of system parameters. There are also review papers 
0]- 0] where various aspects of physics of polarized BEC 
were considered. In this paper we interested in possibili- 



ties of polarization waves existence in DBEC, i.e. another 
type of EE. Dynamics of polarization is interesting not for 
quantum gases only, but in solid state physics and physics 
of low-dimensional systems too. The polarization waves 
in low dimensional and multy-layer systems of conduc- 
tors and dielectrics are considered in the papers I23il24ll . 
Analogously, in BEC of atoms with EDM we expect the 
existence of polarizations wave; along with Bogoliubov's 
mode. 

For investigation of BEC, in system of particles with 
EDM or magnetic moments, various theoretical methods 
are used. One of the ways of theoretical description of the 
polarized BEC is the generalization of well-known Gross- 
Pitaevskii (GP) equation i25ll - lE9Tl . In this way next pa- 
pers were made: for the spinor BEC l3Tj| - [32ll . for the effect 
of magnetic moment on the BEC evolution B2I1 - 14 ill and 
for the influence of electrical polarization of atoms on the 
dynamic processes in BEC [32]- [42]. Dipolar fermions 
also attract a lot of attention fl43|], Hi . 

The same generalization of GP equation usually used for 
magnetic moments and electric dipoles: 
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In equation fl} following designation are used: <£(r, t) is 
the macroscopic wave function, is the chemical potential, 
V ext is the potential of external field, g is the constant of 
short-range interaction, d is the dipole electric moment of 
single atom, m is the mass of particles and H is the Planck 
constant divided by 2tt. 

Along with GP equation the method of quantum hydro- 
dynamics (QHD) is also used. The method of QHD n 1411 . 
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1 4511 . [46] was used for investigations of various physical 
systems, particularly for BEC 11411 . II47I1 and polarization 
waves in conductors and semiconductors l23ll . Therefore, 
the method of QHD could be used to consider the possibil- 
ity of arising polarization waves in BEC and calculating of 
the dispersion law of these waves. 

At interpretation of DDI in quantum gases (QG) used 
the notion of scattering length (SL) and first Born approx- 
imation (FBA), analogously to the SRI described by the 
fourth term in right-side of equation CO. Therefore, the 
problem reduces to the process of scattering. Different ap- 
proximations based on scattering process was analyzed in 
Ref. [48]. Particularly, where were considered condition 
there can be used FBA and presented generalization of GP 
equation for the scattering of polarized atoms beyond FBA. 

There we consider QHD we do not formulate the scat- 
tering problem. We do not need to confine ourselves to a 
scattering process, because we interesting in interaction be- 
tween atoms including interaction between several atoms at 
the same time. This statement is especially important for 
polarized atoms because of the long-range interaction be- 
tween dipoles. Of course for system of particles with low 
concentration each act of interaction between dipoles could 
be considered as a scattering process. This approximation 
is possible as for neutral polarized particles as for charged 
particles. The kinetic theory for the last case was devel- 
oped by L. D. Landau in 1936 14911 . Landau constructed 
collision term for Coulomb particles, using analogy with 
the Boltzmann equation. Landau collision integral depend 
on scattering cross section. In this connection he consid- 
ered the problem of scattering of charged particles on small 
angle, that is suitable for rare systems. Strictly speaking, 
for particles with long-range interaction the conception of 
self-consistent field is more suitable. This conception was 
suggested by A. A. Vlasov in 1938, for system of charged 
particles 15011 . In fact group of particles moves in the field 
of their neighbors, but state of motion of neighbors depend 
on motion of this group of particles as well (This picture 
of interaction is more detailed described on Fig. [[]). Ac- 
cording to this idea Vlasov proposed kinetic equation for 
charged particles called Vlasov's equation. Actually the 
Landau collision integral and the Vlasov's equation are the 
two opposite branches of hypothetical general theory. In 
this paper we follow to Vlasov's approximation and derive 
the self-consistent field theory for dipole-dipole interaction 
in quantum gases. From this point of view there no need to 
consider the first Born approximation or scattering process 
as well. In Appendix B to this paper we present derivation 
of GP equation for polarized particles from QHD equa- 
tions. We did it at the condition that dipoles are parallel 
each other and to external field, they directions does not 
oscillate or vibrate around equilibrium direction. 

Other methods were also applied for investigation of the 
polarized BEC and other systems of particles, with EDM, 
along with GP equation. A hydrodynamic formulation of 
the Hartree-Fock theory for particles with significant EDM 
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FIG. 1. (Color online) The figure presents the picture of self- 
consistent interaction between dipoles. Total polarization of re- 
gion 2 (red) interacts with total polarization of region 1 (green). 
Changing the extreme point of a radius vector (or shifting the 
region 2) we can scan whole space. In this way we obtain ac- 
tion of external dipoles on region 1 . Changing position of region 
1 and repeating described operation we obtain action surround- 
ing dipoles on each region of space. This is a picture of self- 
consistent interaction in fixed moment of time and this picture 
governs an evolution of polarization in system. This picture of 
interaction is typical for classic physics, where we need to obtain 
smooth functions describing collective motion. For that is neces- 
sary to average at physically infinitesimal volume (sketched cir- 
cle). In quantum mechanics, where concentration, polarization, 
etc, defined via wave function and we can consider described pic- 
ture on interaction of separate particle instead of space regions. 



is considered in l5lll . In paper l5lll the Euler-type equa- 
tion was derived, from the evolution of the density ma- 
trix. The EDM dynamics in dimer Mott insulators causes 
the rise of the low -frequency mode [52]. A two-body 
quantum problem for polarized molecules is analyzed in 
115311 . The existence of states with spontaneous interlay er 
coherence has been predicted in |5411 in systems of po- 
lar molecules. Calculations, in 115411 . were based upon the 
secondary quantization approach, where the Hamiltonian 
accounts for the molecular rotation and dipole-dipole in- 
teractions. Superfluidity anisotropy of polarized fermion 
systems was shown in 115 511 and their thermodynamic and 
correlation properties were investigated 156J]. The effect of 
EDM in a system of cooled neutral atoms which are used 
for quantum computing and quantum memory devices was 
analyzed in 115 711 . 

The characteristic property of BEC in a system of ex- 
citons inducing in semiconductors 115 811 it is a significant 
value of EDM of excitons. This leads to strong interac- 
tion of excitons with external electric field and emergence 
of the collective DDI in exciton systems. QHD method 
may be applied to such systems along with quantum kinet- 
ics based on the nonequilibrium Green functions I159ll or 
density matrix equations llrSbTl . 

Notable success has been reached in the Bose condensa- 
tion of dense gases |61, 6^1 . Consequently, we need the 
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FIG. 2. (Color online) The figure shows motion of dipoles with- 
out changing of direction of dipoles during motion. Dipoles stand 
at temperature equal to zero T — and in external uniform static 
electric field E. At lower picture describe particles in two differ- 
ent moments of time. White dipoles presents system in the first 
moment of time when perturb only one dipole, but violet dipoles 
are present system in the next moment. 
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FIG. 3. (Color online) The figure presents the change of polar- 
ization density at movement of dipoles in the absence of dipoles 
direction evolution. 



detailed account of short range interaction in investigations 
of dynamics of EDM. In this work we account the SRI up 
to the third order on interaction radius (TOIR) [14]. This 
approximation leads to nonlocality of SRI B , 1631- 1661. 

Electrically polarized BEC can interact with the beam 
of charged and polarized particles by means of charge- 
dipole and dipole-dipole interaction. Such interaction leads 
to transfer of energy from beam to medium and, conse- 
quently, to generation of waves. In plasma physics the ef- 
fect of generation of waves by electron [67] or magnetized 
neutron [68, 69] beam are well-known. When we use term 
beam we also mean stream of particles excited in the con- 
sidered system, along with external beam of particles pass- 
ing through the system. In presented paper we consider 
similar effect in polarized BEC. 

Our treatment we derived and used the QHD equations 
for polarized particles. System of QHD equations con- 
sists of continuity equation, momentum balance equation, 
equation of polarization evolution and equation of polariza- 
tion current evolution. We made first principles derivation 
of this equation. For this purpose we used many-particle 



Schrodinger equation. In this article we analytically calcu- 
late the dispersion properties of polarized BEC. We have 
shown that the dynamics of EDM leads to existence of 
new branch in dispersion law. Consequently, there are the 
waves of polarization in polarized BEC, along with Bo- 
goliubov's mode. We obtain appropriate contribution to 
the dispersion of Bogoliubov's mode related to DDI. Then, 
we consider the process of wave generation in polarized 
BEC by means of monoenergetic beam of neutral polar- 
ized particles. We predict new method of generation of 
Bogoliubov's mode and a new type of modes, rest of the 
polarization waves. 

This paper is the extension and result of processing of 
our previous paper [70], where the brief description of the 
same results and ideas were presented. 

This paper is organized as follows. We introduce the 
model Hamiltonian in Sec. II, and present the momentum 
balance equation. Farther, in Sec. II, we derive the equa- 
tion of evolution of polarization and obtain the influence of 
the interactions on the evolution of polarization. In Sec. Ill 
we calculate the dispersion dependence of EE in BEC. The 
polarization evolution is taken into account. We obtain the 
contribution of polarization in dispersion of Bogoliubov's 
mode and show the existence of new wave solution. We 
consider high frequency excitation appearing at large equi- 
librium polarizations. We also present calculations for the 
wave of polarization at constant concentration. In Sec. IV 
we study the wave generation in polarized BEC by means 
of neutral polarized particle beam. In Sec. V we present 
the brief summary of our results. App. A. contains de- 
tail of derivation of equations of polarization evolution. In 
App. B. derivation of GP equation for polarized particles is 
presented. In App. C. we consider distinction of evolution 
of spinning particle from the particles with EDM. 



II. THE MODEL 

In GP equation term described DDI, for both EDM and 
magnetic moment, interpret via dipole scattering length 
(DSL) add- The angle between dipoles coming in GP equa- 
tion via potential of DDI, and equation is solved for various 
angles between dipoles. But equation does not contain in- 
formation about evolution of dipoles direction or, for the 
case of soft dipoles than dipoles caused by external field, 
variation of dipole magnitude. Change of dipoles direction 
there is due to interaction. If we have external static uni- 
form electric field, at close to zero temperatures, all dipoles 
will be direct parallel to external field. At low tempera- 
tures there is no temperature disorientation of dipoles and 
we have system of parallel dipoles, as it is shown on Fig. 
121 In such systems there are two type of dipoles motion. 
The first one is the motion there the direction of dipoles 
is preserves Fig. |2 [3] In this case DDI has influence on 
dipoles properties because it is change the strength of in- 
teraction and leads to anisotropy of interaction. At such 
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FIG. 4. (Color online) The figure describes how the polarization 
direction perturbation propagates through the system of particles. 
We consider system of dipoles at T = and being in external uni- 
form electric field. Upper picture show that if at the first moment 
there is perturbation of the EDM in next moment this perturbation 
pass to several neighbors. At lower picture describe particles in 
two different moments of time. White dipoles presents system in 
the first moment of time when perturb only one dipole, but violet 
dipoles are present system in the next moment. 



form of motion the polarization is changed accordingly to 
concentration P(r, t) = dn(r, t), where d is the EDM of 
single particle Fig. [3] The second type of motion is a mo- 
tion of dipoles including change of dipoles direction Fig.s 
01 |5j Perturbation of dipoles direction can be caused by 
both external influence and local fluctuation of dipoles di- 
rection produced by DDI. In consequence of DDI the per- 
turbation will propagate through system. Particularly, this 
propagation might be the wave nature. For description of 
such motion we need to know how polarization of system 
evolves in different point of space during the time. There- 
fore we must derive equation describing the evolution of 
polarization field and accounting influence of interaction 
of this evolution. 

Angle dependence in GP equation is important for the 
case when external electric field is nonuniform. In this case 
dipoles in different areas of space have different directions 
Fig. [6] and it is reflected on interaction via angle depen- 




FIG. 5. (Color online) The figure presents the oscillation of EDM 
near of the equilibrium position which defined by external elec- 
tric field E. Indexes 1-6 present the different moments of time 
ti < <a < •■• < tg. This figure describes mechanism of the 
polarization change due to oscillations of EDM of each particle 
in the system. This mechanism might reveal along with rotation 
of EDM. Incoherent motion of EDMs which is combination of 
rotation and oscillation is presented on Fig. |5b). 



FIG. 6. (Color online) The figure presents the two dipoles with 
the EDM turned in different directions. At propagation of per- 
turbation through the system relative directions of dipoles is 
changed. This lied to change of force of interaction in concor- 
dance with potential of DDI (O. In Appendix B we shown that 
GP equation describe motion for parallel dipoles only. 



dence under integral in GP equation £Q|. Still, this equation 
does not contain information about propagation of dipoles 
direction perturbation. 

For studying we pay attention for the SRI and DDI. We 
also consider the action of external fields on particles. The 
action of external electric field on dipoles we consider in 
evident form and the external field described trap presented 
in general form by trapping potential V tra p(^i, t). 

Here we used quasi-static approximation for DDI. The 
meaning of quasi-static approximation is that we consider 
interaction between dipoles at fixed moment of time as they 
are motionless. DDI leads to change of state of motion, 
and, in particularly, to variation of dipoles direction. Then 
we talk about change direction we mean both the varia- 
tion during precession around of rotation axis and devia- 
tion of rotation axis direction. Described approximation is 
an analog of Coulomb interaction there interaction between 
moving charged approximately described by means of in- 
teraction potential of motionless charges, this approxima- 
tion usually used for description of non-relativistic plasma 
& 

Explicit form of the Hamiltonian of considering system 
in a quasi- static approximation is 
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2m. 



-2 



(2) 



The first term here is the operator for kinetic energy. The 
second term represents the interaction between the dipole 
moment df and the external electrical field. The subse- 
quent terms represent short-range and dipole-dipole 
interactions between particles, respectively. The Green's 
function for dipole-dipole interaction is taken as Gff = 

V?Vf(l/ry). 

For investigation of EE dynamic in polarized BEC we 
derive the system of QHD equations. This system of equa- 
tions consists of the continuity equation, the Euler's equa- 
tion and for the case of polarized particles the equations of 
polarization evolution and equation of field. The system of 
equations is derived by methods described in ifbUl . 

The first equation of a QHD equations system is the con- 
tinuity equation 



d t n(r, t) + d a (n(r, t)v a (r, t)) = 0. 



(3) 



The momentum balance equation for the polarized BEC 
has the form 

mn(r, t)(d t + vV)/(r, t) + d p af3 (r, t) 

-JLr Mr, t ) + ( ^4^M) 

Am 4m \ n(r,t) J 



Tn(r,t)d a n(r,t) + -T 2 d a An 2 (r,t) 
+ P p (r,t)d a E l3 (r,t), 



where 



and 
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Or 



dr{r) 
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(4) 



(5) 
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30 J ' v ' dr 

In equation (© we defined a parameter T 2 as This 
definition differs from the one in the work ifbUl . Terms 
proportional to Jl 2 appear as a result of usage of quantum 
kinematics. The first two members at the right side of the 
equation (O are first terms of expansion of the quantum 
stress tensor. They occur because of taking into account of 
the SRI potential Utj . The interaction potential Uij deter- 
minates the macroscopic interaction constants T and T 2 . 
The last two members of the equation Q describe force 
fields that affect the dipole moment in a unit of volume as 



the effect of the external electrical field and the field pro- 
duced by other dipoles, respectively. The last member is 
written using the self-consistent field approximation II45I1 . 
l23ll . p a P(r,t) is a tensor of the kinetic pressure, which 
depends on particles' thermal velocities and does not con- 
tribute into the BEC dynamics at temperatures near zero. 

The first order of the interaction radius interaction con- 
stant for dilute gases has the form 



T 



4irh 2 a 



in 



(7) 



where a is the scattering length SL HQ. The value T 2 
may be expressed approximately as 

0a 2 T 

T 2 = — , (8) 

where is a constant positive value about 1 , which depends 
on the interatomic interaction potential. Finally, T 2 takes 
the form 

ir9h 2 a 3 

T 2 = . (9) 

2m 



(10) 



We have also a field equation 

VE(r,t) = -4vrVP(r,t). 

In the case particles does not contain the dipole mo- 
ment, the continuity equation and the momentum balance 
equation form a closed system of equations. When the 
dipole moment is taken into account in a momentum bal- 
ance equation, a new physical value emerges, a polariza- 
tion vector field P a (r, t). This causes system of equations 
to become incomplete. 

Next equation we need for investigation of EE dispersion 
is the equation of polarization evolution 

d t P a {r,t) + d p R ap (r,t) = 0, (11) 

R a P(r, t) is the current of polarization. 

The equation (fTTT l does not contain information about 
the effect of the interaction on the polarization evolution. 
The evolution equation of R a/3 (r, t) can be constructed by 
analogy with the above derived equations. Method of the 
equations derivation is described in Appendix A. Using a 
self-consistent field approximation of the dipole-dipole in- 
teraction we obtain an equation for the polarization current 
R a0 (r,t) evolution 

d t R aP (r, t) + d 1 [R al3 (r, ty (r, t) 



+R a ~<(r,t)v p (r,t) - P a (r,t)vP(r,t)v^(r,t) 



m Am z 

, h2 dl ( dpP a (r,t)d in (r,t) | d 7 P"(r, t)d p n( r , t) 
8m 2 \ n(r,t) ' n(r,t) 
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— Td p (n(r,t)P a (r,t)] 
m \ J 



+ i ,P»(M)^(M) 
m n(r,t) 



(12) 



Here r a/37 (r, i) represents the contribution of thermal 
movement of polarized particles into the dynamics of 
R al3 (r,t). As we deal with BEC below, the contribution 
of r Q/37 (r, t) may be neglected. The last term of the for- 
mula ([T21 includes both external electrical field and a self- 
consistent field that particle dipoles create. This term con- 
tain numerical constant a. Used in formula (fT2l approxi- 
mations are described in Appendix A. 

The first term in right side of Eq. (fT2l) describe the short 
range interaction. 

We can see various interactions are included in the equa- 
tions © and dTzb additively. At a short distances among 
particles acts both SRI and dipole-dipole interaction. At 
large distances remain only dipole-dipole interaction. 

From Eq. (fT2l we can see that the change of polariza- 
tion arise from both the dipole-dipole interaction and the 
short range interaction. The SRI among particle leads to 
displacement of particles. Consequently, as particles has 
EDM, where are motion of EDM, e.i. changing of the 
i? Q/3 (r,t). 

The terms that are proportional to h 2 are of quantum ori- 
gin as they are analogs of the Bohm quantum potential in 
the momentum balance equilibrium. 

In the Appendix C we consider the difference between 
equations of electrical and magnetic polarization evolution. 

QHD of spinning particles with spin-orbit interaction 
were developed in Ref. l69ll for system of charged parti- 
cles, however where obtained field of the force can be used 
for neutral quantum gases too. 



IU. ELEMENTARY EXCITATIONS IN THE POLARIZED 

BEC 

We can analyze the linear dynamics of elemental exci- 
tations in the polarized three dimensional BEC using the 
QHD equations ©, ©, (QHl, CO and ([Hi. Let's as- 
sume the system is placed in an external electrical field 
E = E e z . The values of concentration n and polar- 
ization P = kE for the system in an equilibrium state 
are constant and uniform and its velocity field v a (r, t) and 
tensor R a/3 (r, t) values are zero. 

We consider the small perturbation of equilibrium state 
like 



n 



n + 5n, v a = 



p a = pa + 6 pa Ra _ _|_ $ Ra p_ 



(13) 



Substituting these relations into system of equations (O, 
©, (TTTT ). ( fT2l) and dTOl) and neglecting nonlinear terms, we 



obtain a system of linear homogeneous equations in par- 
tial derivatives with constant coefficients. Passing to the 
following representation for small perturbations 5f 

5f = f(co, k)exp(—iujt + ikr) 

yields the homogeneous system of algebraic equations. 
The electric field strength is assumed to have a nonzero 
value. Expressing all the quantities entering the system 
of equations in terms of the electric field, we come to the 
equation 

A • E K = 0, 



where 



A = co 2 



h 2 k 4 Tk 2 n . P 2 k 2 



4m' 



+ 



2m 



+ 4tt<7-^ 



mn 



27rTk i P 2 



m 2 ui 2 — h 2 k 4 /4 + mTk 2 n — mT 2 k 4 n 
In this case, the dispersion equation is 
A = 0. 

Solving this equation with respect to uJ 2 we obtain a fol- 
lowing results. 

The dispersion characteristic for EE in BEC can be ex- 
pressed in the form of 



1/3 h 2 k A 
oj 2 = —\ --Tn k 2 + -— + T 2 n n k 4 + 4trt 

2m \ 2 2m n Q 



P 2 k 2 




Tn Q k 2 — T 2 n k i + 47r<7- 



P 2 k 2 



8irTk 4 P 2 



(14) 



In contrast to the non-polarized BEC, where only Bo- 
goliubov's mode exists 0, B, a new wave solution 
appears in a polarized system due to the polarization dy- 
namics. Bogoliubov's mode corresponds to a solution of 
the equation (TBI with minus sign before the square root. 
New wave solution it is a wave of polarization. In general 
case presented by formula (TBI the frequency of polariza- 
tion wave in BEC depend on P and T. To investigate 
the BEC polarization effect on the Bogoliubov's mode and 
the dispersion characteristic of the new solution we analyze 
extreme cases of the formula (fT4ll . 

Formula ([LUl demonstrates that taking account of the 
BEC polarization dynamics leads to a new solution. 

Let's start with the case when contributions of different 
term in formula (fT4l l are comparable. In this case we stud- 
ied dispersion numerically. For this purpose we represent 
formula (fl4ll in the following form 



e 



k 2 



3 a 

4 \a\ 167rn |a| 



da 3 k 2 6vrcra 



16 a 
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1 a 9a 3 k 2 Giraadd 
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FIG. 7. (Color online) Here we use lg k, where k measured in 
cm -1 . The figure presents the dependence reduced frequency £ 
on the wave vector k and DSL <z<2d/M f° r Bogoliubov's mode. 
Figure was obtained for repulsive SRI a > 0, T < 0, \a\ 



10" 



cm, no = 10 cm 




FIG. 8. (Color online) Here we use lg k, where k measured in 
cm -1 . The figure presents the dependence reduced frequency £ 
on the wave vector k for polarization wave at three different value 
of DSL a dd : a ddA = (red), a dda = \a\ (green), a dd>3 = 10|o| 
(blue), a dd 4 — 100|a| (black). Figure show increasing of £ with 
increasing of a dd . Figure obtained for repulsive SRI a > 0, T < 
0, and \a\ = 10~ 7 cm, n = 10 14 cm~ 3 , a = 1, 9 = 1. 



where 



e 



\T\nok 2 /m 

Here we use the length a dd measures the strength of dipolar 
interaction, and defined by formula 



a dd 



Yl'Knlh 2 ' 



Above we discuss that we do not consider the scattering 
process, but for convenience and analogy with the results 
of other paper we enter this parameter, although a dd does 
not connect with scattering. 
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FIG. 9. (Color online) Here we use lg k, where k measured in 
cm" 1 . The figure shows behavior of reduced frequency £ depen- 
dence on wave vector k for both solutions presented by formula 
d 1 5 b - We can see the relative behavior of two solutions. Upper 
curve (blue) coincides to solution with sign "+" in front of square 
root (polarization mode) and lower curve (red) correspond to the 
one with sign "-" (Bogoliubov's mode). a dd = \a\, a > 0, T < 
(repulsion), \a\ — 10~ 7 cm, no = 10 14 cm -3 , <r = 1, 9 = 1. 
The figure presents comparison of the reduced frequencies de- 
pendence for Bogoliubov's and polarization modes. We can see 
that £(fc) for polarization mode lie above Bogoliubov's mode. 
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FIG. 10. (Color online) Here we use lgfc, where k measured in 
cm" 1 . The figure presents the dependence reduced frequency £ 
on the wave vector k for Bogoliubov's mode, a < 0, T > 
(attraction), DSL a dd = (Red, the right one), a dd = 0.0005|a| 
(Blue), a dd = 0.005|a| (Green), a dd = 100|a| (Red, the left one), 
\a\ = 10~ 7 cm, n = 10 14 cm -3 , a = 1,6 = 1. 



Here we describe the numerical analysis of formulas 
( fT5l) . Reduced frequency £ for two waves presented by 
formula dT5l) exhibit on Fig.s|7]and[8j They are obtained 
for repulsive SRI (a > 0, T < 0) at a = l(T 7 cm 
and n = 10 14 cm -3 . These Fig.s show the reduced fre- 
quency for different values of equilibrium polariza- 
tion and, consequently, for the different DSL a dd . Particu- 
larly, where the limit case a dd is considered too. 

Fig- E] present solution with sign "-" in front of square 
root in formula (TBI) , this solution corresponds to the Bo- 
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FIG. 11. (Color online) Here we use lgfc, where k measured in 
cm -1 . The figure presents the dependence reduced frequency £ 
on the wave vector fc and DSL a dd for polarization mode, a < 0, 



T > (attraction), 
6 = 1. 



\a\ = 10 cm, no 



= 10 14 cm -3 , 



(7 = 1, 
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FIG. 12. (Color online) Here we use lg k, where k measured in 
cm -1 . The figure presents the dependence reduced frequency 
£ on the wave vector k for polarization mode, a < 0, T > 
(attraction), DSL add = (Red), add = \a\ (Blue), add — 10|aJ 
(Green), a dd = W0\a\ (Brown), \a\ = 10~ 7 cm, n = 10 14 
cm' 3 , a = 1,0 = 1. 

goliubov's mode. For all a d d/\a\ the do not depend 

on a dd /\a\. 

Fig. [8] show strong dependence of reduced frequency 
of polarization mode ("+" in front of square root in ( TBI) ) 
on equilibrium polarization. On FigJS] as on other Fig.s 
this dependence present via dependence on DSL divided 
by module of SL of SRI a d d/\a\- 

The form of £(k) for two solutions has similar form. We 
present both solutions on Fig. [9] at fixed add/\o\ = 1 to 
compare. From Fig. |9]we see that reduced frequency of po- 
larization mode is larger than one's for Bogoliubov's mode. 

From the Figj8]we obtained for polarization mode that 
polarization lead to growing of the reduced frequency £. At 
small wave vectors of order 10 4 cm _1 and for the add = \o\ 
the reduced frequency £ increases in 6 times. For the 
o-dd = 100|a| one's increases in 60 times. With increas- 



FIG. 13. (Color online) Here we use lgfc, where fc measured in 
cm" 1 . The figure presents the dependence reduced frequency £ 
on the wave vector fc for polarization mode. On this Fig. the 
region of instabilities is also shown. It is a region where a gap in 
curves. Figure obtained for attractive SRI (a < 0, T > 0), a d d = 
(Red), add = 0.1|a| (Blue), add = 0.2|a| (Green), add = M 
(Black), \a\= 10~ 5 cm, n Q = 10 14 cm -3 , a = 1,6 = 1. 



ing of wave vector the influence of polarization became 
smaller. For example at k = 10 6 cm _1 and add = 100|a| 
the reduced frequency £ increase just in 1.6 times. 

For the case of attractive SRI in the absence of polar- 
ization there is instability of Bogoliubov's mode at small 
wave vectors. At the presence of small polarization the 
Bogoliubov's mode remains unstable, but area of stability 
become some wider as it shown on Fi sfTUl There is fast 
widening of a region of stability at increasing of a d d/\o\ 
FigHU Widening starts at approximately add = 0.02|a|. 

For the case of attractive SRI in the absence of polariza- 
tion there is instability of polarization mode at small wave 
vector's (see Fig. fT2l >. At a d d = we see instability of 
waves at wave vector k below 1.6 • 10 4 cm -1 . Existence 
of EDM leads to increasing of £ and fast growing of the 
region of stability. Wave are stable up to A; = 10 3 cm -1 
at a-dd/ 1 o\ = 1- Further increasing of a d d/\a\ leads to in- 
creasing of £ up to 64. £ become equal to 64 at add/\a\ — 
100 and does not change at further increasing of add/\o\. 

Comparing Fig.s [TOl and [T2l we can see that reduced fre- 
quency £ of polarization mode much larger than Bogoli- 
ubov's mode, as it was for repulsive SRI. 

With increasing of strength of SRI the role of T 2 become 
more important. It's influence reveals at \ a\ = 10~ 5 cm for 
attractive SRI. In this case there is instability region at rel- 
atively large wave vectors. For described case the reduced 
frequency is presented on Fig.s [L3l and [T4l at various 
&dd/\a\- The region of instability is shown where. This 
instability appears because the expression under the square 
root in formula (TBI) became negative as it shown on Fig. 

El 

Below, we analytically consider limit cases of general 
dispersion solution (fT4b . 

Here we consider the case when the effect of polarization 
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FIG. 14. (Color online) Here we use lg k, where k measured in 
cm -1 . The figure presents the dependence reduced frequency £ 
on the wave vector k for Bogoliubov's mode, a < 0, T > 
(attraction), DSL add = (Red), add = 0.1|a| (Black), add = 
0.2|o| (Blue), a dd = 2|«x| (Green), \a\= 1(T 5 cm, n = 10 14 



Here we use indexes "B" for Bogoliubov's mode and "P" 
for new polarization mode. When deriving formulae ([TBI 
and ( fTTl l we expanded a sub-radical expression in (fT4l and 
took only first two terms of the expansion to estimate the 
influence of polarization on the wave dispersion. 

As it follows from equations (fTTTl and Sl2\ . changes in 
polarization can occur due to dipole-dipole interactions, 
SRI and quantum Bohm potential, i.e. to members propor- 
tional to ti 2 . That's the reason for existence of other ways 
of polarization change when the contribution of equilib- 
rium polarization P is negligible. In the latter case, taken 
in a linear approximation, changes in the polarization do 
not affect the concentration evolution. 

Formulae (TBI) are valid even in the absence of the ex- 
ternal electrical field when equilibrium polarization equals 
zero P = 0. Equations (fT4l in that case take the form 



h 2 k 4 

Am 



uj r = — I 



Tn Q k 2 + T 2 n Q k 4 



(18) 
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FTG. 15. (Color online) The figure presents the behavior of the 
square root in formula d31| > (we conditionally designated this 
quantity as "sqrr") at the change of the wave vector k correspond- 
ing to the Fig. [T0]at a d d = 0.1|a|, (square does not depend on 
concentration). |a|= 10 -5 cm, a = 1, 6 = 1. Here we use Igk, 
where k measured in cm -1 . 



is low compared to the contribution of short-range effects, 
i.e. when the members proportional to T and Y 2 is com- 
parably large. If so, then formula (fT4] l takes the form 

, h 2 k 4 



Am 2 



Tn k 2 ^ T 2 n k 4 



in 



in 



AnP 2 k 2 T(cj - 1) - 2aT 2 k 2 



and 



mn 



h 2 k 4 

Am 2 



T - 2T 2 k 2 

1 Tn fc 2 

2 m 



(16) 



+ 



AirP 2 k 2 T(<7 - 1) - 2oT 2 k 2 



mn Q 



T 



(17) 



U), 




Tn k 2 



(19) 



The waves of polarization could be existed at the absence 
of external electric field E = 0. In this case the equilib- 
rium state is no polarized and dipole direction of particles 
is distributed accidentally. 

If the contribution of equilibrium polarization into BEC 
dispersion (fT4l l is comparable to the contribution of SRI 
in the third order of the interaction radius T 2 relationships 
(fT4l) are transformed into 



h 2 k 4 Tn k 2 T 2 n k 



Am 2 



+ 



+ An- 



ni 



rn 



p2p 

mn 



and 



h 2 k 4 1 Tn k 2 



Am 2 



+ Air- 



rn 



P 2 k 2 
mn 



(a-1) 



(20) 



(21) 



Using Feshbach's resonance 11711 172H we can transform 
the SRI potential in such a way that T = while T 2 ^ 0. 
Formulas (TBI) in this situation turn into 



ui. 



h 2 k 4 T 2 n fc 4 
Am 2 m 



and 



21,4 



K z k 
Am 2 



2 1.2 



+ Aira : 



mn 



(22) 



(23) 



In this paper we primarily focus on the influence of BEC 
polarization on its dispersion characteristics. So, let's con- 
sider the case when the contribution into the dispersion of 
the SRI at the first order of the interaction radius, i.e. terms 
proportional to T, is comparable to the contribution of po- 
larization, and their total effect is much greater than the 
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contribution of terms proportional to T 2 . Here, (fT4l) turns 
into 



and 



h 2 k 4 Tn a k 2 /(a-l 



4m 2 



m 



a 



, H 2 k 4 A P 2 k 2 
4m z mn 

1 Tn k 2 {{(7 + 4) 



■m 



2a 



(24) 



(25) 



IV. GENERATION OF WAVES IN POLARIZED BEC 

In that section we consider the process of wave genera- 
tion in BEC by means of the beam of neutral polarized par- 
ticles. The interaction between beam and BEC has dipole- 
dipole origin. 

To get the dispersion solution we use the system of QHD 
equations for each sort of particles J3J, (|4), (fTTT l. (fT2l and 
the equation of field (flOl) . The equilibrium state of system 
is characterized by following values of the BEC parame- 
ters: 



n 



n + Sn, v a = + v a , 



polarization evolution at constant concentration 

DBEC usually is the rare gas and we can not use condi- 
tion incompressiblity, but we consider the case there con- 
centration of particles is not change in time. We interested 
in dynamics of wave of polarization only. For this aim we 
need to use equations (fTTT l and (fT2| at additional condition 
n(r, t) = n = const. This is means the polarization 
change due to evolution of dipoles direction. Equation of 
polarization evolution (fTTT) has no change in this case. In 
eq. ( fT2l) nonlinear term are disappeared and the term de- 
scribing the SRI is simplified: 

d t R ap (r, t) + d"> ( R al3 (r, t)v~<(r, t) 



pa = pa + § p a i R a0 = Q + SRa yg) 

and values of the beam parameters: 

n b = n ob + 5n b , v% = U5 za + 8v£, 



TDOt 



P a b + SP b a , R 



a/3 



R^ + SRf. 



(29) 



The polarization Pq is proportional to external elec- 
tric field E£. We consider the case then E = 
[E sin<p, 0, E cos<f]. In this case the tensor R^ has 



only two unequal to zero elements: R, 



zx 
Ob 



R ob simp and 



i?ob = Ro b costp. For the process, under consideration the 
dispersion relation is: 



1 + 4irk 2 i 



p2 
r 



ft 2 fc 4 , Tn k 2 
Am 2 ' 2m 



+iT 7 (r, t)v (r, t) - P a (r, t)v p (r, t)^(r, t) 



d p AP a {r,t) = ^Td^f P a (r,t)) 
m \ J 



4m 2 



Tk 2 /{2m 2 ) 



a 



h 2 k 4 i Trtpfc 2 
4m 2 1 m 



T 2 n fc 4 



k z u) 2 - m 



mn 



+ —P a (r, t)P 7 (r, tWEnlr, t). (26) 
mn 

and equation of field dTOl) . 

In this approximation exist only one wave solution ex- 
ists, it is the wave of polarization. It's dispersion is 



to 



h 2 k A (4-koP 2 
4m 2 V 



mn 

In this case formula (1271) reduces to 



4iraP 2 



m I 



(27) 



mn 



For the particles with the EMD equal to 1 Debye at con- 
centration 10 12 cm -3 and wave vector k ~ 10 7 cm -1 the 
frequency u) reach value of order 10 6 — 10 7 s _1 . 

Dispersion of ID and 2D wave of polarization is de- 
scribed in Ref. II23II . 



2(u- k z U)P Qb k z (P ob U - R ob ) a b P 2 b \ 
m b n ob \ (w - k z U) 2 - I 



0. 



Using relation P ob U — R Qb 
tion (l30l and obtain 



(30) 

0, we can simplify the equa- 



te 



lo 2 d ( Tn k 2 /(2m) 



a 



°bUJ 2 Db 



ft 2 fc 4 

4m? 



0. 



(w - k z U) 2 

In this formula the following designations are used 

h 2 k A Tn k 2 



4m 2 



2m 



(31) 



(32) 
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2 1.4 



and 



4m 2 



UJ, 



m 



4ttP^ 

m i n Oi 



m 



(33) 



(34) 



where i is the index of sorts of particles, the BEC or the 
beam. 

The equation d3~TT) has two beam related solutions, in the 
absence of BEC medium: 



UJ 



k z U± 



>h 2 k A 

4m\ 



+ °bU 2 Db . 



(35) 



We will consider the possibilities of instabilities for the 
case of low-density beam, the limit case ui Db ~ n ob — > 0. 
In this case we can neglect the last term in square root in 
d35l) . The resonance interaction beam with the BEC realiz- 
ing if 



k z U± 



hk 2 

2m b 



u;(k), 



(36) 



and could lead to instabilities. The quantity ui{k) is the 
dispersion of BEC modes (TPTl . The frequency in this case 
can be presented in the form 



uj = k z U ± h qui. 

2m b 

Let us to consider two limit cases. 

small frequency shift limit 

In the limit case 

5u) <C hk 2 /m 
the frequency shift obtain in the form: 



(37) 



5u 2 



where 



2a b m b m 2 ui 2 Db (ijj 2 — uj 2 ) 2 {uj 2 — wf) 2 
: LJulhr 2 nlk 6 W 



2ma 



(38) 



(39) 



W = 2u? - w? - u 2 2 - y^^ 2 ~ w 2 2 ) 2 - (40) 

and the frequency uj determined with formula (TT4l . The 
instabilities take place in the case 5ui 2 < 0. The sign of 
5u) 2 is depend on the sign of W. 

For the case resonance interaction of beam with the 
waves in BEC there are instabilities, for the first beam 
mode in d35l) at W < and for the second beam mode 
in formula (l35l l if W > 0. For the polarization mode W is 
positive. It means that the interaction of polarization mode 
with the second beam related mode results in the instabil- 
ity. For the Bogoliubov's mode the sign of W depend on 
a. 

We can consider the following cases: 



(i) the contribution of equilibrium polarization to ui{k) 
is dominant; then W > 0; 

(ii) the dominant contribution to ui{k) results from the 
term in ui(k) which is proportional to T, and, equilibrium 
polarization and SPI in TOIR give comparable contribution 
to uj(k). In this case there is 



00 



1 + 



T 2 n k 4 jm 



Tn k 2 /(2m) + 2T 2 n k 4 /m - Aooj 



. The sign of W varies at a = cr . The dependence of W 
from a is presented in the Table 1 . 



large frequency shift limit 



In the limit case 



we have 



5u) 3> hk 2 /m 



(41) 



r t 3/ (r b m 2 uj 2 Db {uj 2 - uj 2 ) 2 (uj 2 - UJ 2 ) 2 
y uiuj z D T 2 nik 4 I W I 

where £ equal to £1 = \/l for W > or = for 
W < 0. Evident forms of quantities £1 and £_i are 



[-v 



-1 + n/3 -1 -iVs, 



and 



6 = [1, 



The considerations concerning the sign of W, reflected 
in the table 1, are valid also for the limit condition d4"TI) . 

From the formulas d39l and d42l) we can see that external 
neutral particles beam leads to instabilities for both Bogoli- 
ubovs and polarization waves. 



TABLE I. In this table the sign of the W are presented for the 
Bogoliubov's mode when SRI is prevailed. 





T>0 


T<0 




attraction 


repulsion 


a > a 


+ 




a < a 




+ 
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V. CONCLUSION 

In this work we developed the self-consistent method for 
description of dynamics of polarized BEC. This method 
accounts for the effect of polarization on changes in the 
concentration and in velocity field, which are determined 
in general by the continuity equation and Euler's equation. 
We derived the evolution equations of the polarization and 
the polarization current. The derived equations contain in- 
formation about the influence of the interactions on the po- 
larization evolution. We studied the effect of polarization 
on the BEC dynamics and the influence of SRI on the po- 
larization evolution. An expression of SRI contribution in 
the equation of polarization current evolution via concen- 
tration, polarization and the SRI potential T = Y(£/y) 
was derived. With the assumption that the state of polar- 
ized Bose particles in the form of BEC can be described 
with some single-particle wave function. Changes in polar- 
ization due to SRI are shown to be determined at the first 
order of the interaction radius by the same interaction con- 
stant that occurs in Euler's equation and Gross-Pitaevskii 
equation. 

The derivation of GP equation for polarized particles 
from QHD equation obtained in the paper. The condition 
of validity of GP equation was presented. Comparison of 
evolution equation of electrically polarized BEC and mag- 
netized BEC was discussed. The differences between these 
equations are shown. 

The physical meaning of self-consistent approximation 
was described. Suitability of self-consistent approximation 
for electrically polarized BEC was shown. Distinction of 
developed and used here self-consistent approximation for 
dipole-dipole interaction from scattering process are dis- 
cussed. 

Correct form of the dipole-dipole interaction Hamilto- 
nian is discussed. The arguments for choosing of correct 
form of Hamiltonian for electric and magnetic dipoles are 
presented. Corresponding correct Hamiltonians were de- 
rived. 

The dispersion of EE in the polarized BEC was analyzed. 
In the article, we show that polarization evolution in BEC 
causes a novel type of waves in BEC. The effect of polar- 
ization on the dispersion of the Bogoliubov's mode and the 
dispersion of a new wave mode were studied. 

We show the possibility of wave generation in polarized 
BEC by means of monoenergetic beam of neutral polarized 
particles. 



APPENDIX A 

Hamiltonian's correct form 

At derivation of system of QHD equations we need to 
write evident form of Hamiltonian of dipole-dipole in- 



teraction. In some works the dynamics of the magnetic 
dipole moment and of the EDM ll32l p74ll are analyzed 
in similar ways. Usual expressions of a Hamiltonian for 
dipole-dipole interaction are equal for electric and mag- 
netic dipoles: 

H dd = 6 - *LLlLdZ4 l (43) 

In this paper for interaction of electric dipoles we used fol- 
lowing Hamiltonian: 



H. 



dd 



r 



There is well-known identity 

_ gagpl = i 6rrjr_ + 4^ 

r r A 6 

so we can see the difference between usually using Hamil- 
tonian and one's used in this paper. The corrections of our 
selection followed from the fact that the equations obtained 
in the paper coincide to the Maxwell equations. At deriva- 
tion of QHD equation in self-consistent approximation lead 
to field equation ([Tot , we rewrite it here 

VE(r,t) = -4vrVP(r,t), (45) 

but if we used Hamiltonian ( 1431 we would obtain 

VE(r,t) = yVP(r,t), (46) 

instead of (l45i 

It has been shown by Breit IT76Tl that a Hamiltonian for 
spin-spin interaction, and, as a consequence, for the inter- 
action of magnetic moments contains a term that is propor- 
tional to Dirac 5-function S(r\ — r 2 )d°d2 along with (l43"T l. 
The coefficient of the <5-function has been refined later I177ll 
so that the Hamiltonian is in accord with Maxwell's free 
equations, such as divB = 0. The resultant expression for 
the spin-spin interaction Hamiltonian is: 



47r$ a/J $(r 12 ) + V?V?(l/r 12 ) U?^, (47) 



<*,3 



or 



Actually, the 5 functional term is used in literature I175ll . 
but authors followed to Breit and Landau [76], they used 
coefficient — 4tt/3 at S function instead of 87r/3, but this 
form of coefficient does not corresponds to Maxwell equa- 
tion. It was shown in Ref. [77], where also were obtain 
correct coefficient presented in formula ( PTTT l. 

Thus, the differences in the dipole-dipole interactions of 
electric dipoles, and magnetic dipoles, must be taken into 
account in the development of theoretical field apparatus. 
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Method of equations derivation 

The Schrodinger equation defines wave function in a 
3N-dimensional configuration space. Physical processes 
in systems that involve large number of bodies occur in 
a three-dimensional physical space [78]. That's why a 
problem evolves of obtaining a quantum-mechanical de- 
scription of a system of particles in terms of material fields 
e.g. concentration, momentum density, energy density and 
other fields of various tensor dimension that are defined in 
a three-dimension space. 

The first step in MPQHD equation derivation is the defi- 
nition of particles concentration. We define the quantum 
particles concentration as a quantum average of classic 
microscopic concentration in coordinate representation on 
many-particle wave function ijj(R,t). Thus, the particles 
concentration evident form to be 

n(r,t)= [ dRj2 5(r-ri)ip*{R,t)ip(R,t), 



where dR = YliLi ^ r »- Differentiating this function over 
time and using the Schrodinger equation with the Hamil- 
tonian 10, we derive the continuity equation ((3). In the 
continuity equation the new quantity appears. We can dif- 
ferentiate this quantity over time and use the Schrodinger 
equation with the Hamiltonian (O. In the results we ob- 
tain equation for this quantity. In this way we can find an 
evolution equation for every additive physical quantity, but 
we mast know it's definition. The definition of new phys- 
ical quantities appears at derivation of QHD equation. For 
example, the probability current definition appears at the 
derivation of continuity equation. At derivation of the cur- 
rent evolution equation (O appears several quantities, one 
of them is the polarization presented below (l48l l. Knowl- 
edge of the polarization definition give us ability to derive 
the equation of polarization evolution ( PT2l) and so on. 
In equations (©, ( fTOt polarization occurs in the form of 

P a {r,t) = [ dRJ2S(r-*i)r(R,t)dW(R,t), (48) 

J i 

where r< is the coordinate operator for i-th particle, dR = 
IlJLi dr P- 

In right hand side of equation ( PT2l) where is the force-like 
field F a P(r, t) which give rise to evolution of the polariza- 
tion current R afi (r, t). In general form, for F a @(r, t), we 
can write: 

F a/3 {r,t) = -—d^ a ^(r,t) +—D a Vr,t)d p E' l (r,t). 
m m 

(49) 

A tensor 

D a0 (r,t)= [ dRj^Hr -r t )dy^*(R,t)ij(R,t), 

(50) 

occurs in the term that represents a dipole-dipole interac- 
tion and an interaction of the dipole with external electrical 



field. Based on the reasons of dimensions this value can be 
approximately presented as 

D . , )= « 

n(r,t) 

A SRI causes the tensor E a/37 (r, t) to occur in the equa- 
tion (fl2l> . Taken at the first order of the interaction radius it 
has the form 

E^(v,t) = -^JdRj2S(r-R ij ) 



f^.fl. flJJ. . 



r 



i j 



dr 



(52) 



Tensor £ Q,37 (r, t) describe the influence of SRI on evolu- 
tion of polarization. Eq. (l52l describe the SRI. 



If we apply the procedure described in 11411 to the cal- 
culation of the quantum stress tensor, and neglect the con- 
tribution of thermal excitations, £ Q ^ 7 (r, t) for BEC takes 
a form of 



BEC 



M) 



i 



T5^n{r,t)P a (r,t). (53) 



Formula (1531 is obtained for the case where particles lo- 
cated in state with the lowest energy, which could be de- 
scribed by one particle wave function. This state may be 
the product of strong interaction. 

Tensor £ Q/37 (r, t) is, therefore, like the quantum stress 
tensor a al3 (r,t) in the momentum balance equation ©, 
dependent on T at the first order of the interaction radius 

iH. 



APPENDIX B: DERIVATION OF THE NON-LINEAR 
SCHRODINGER EQUATION 

In Ref. ifbUl the derivation of GP equation was per- 
formed for non-polarized particles from QHD equation, 
here we present one for polarized particles. The NLSE 
comes about from the continuity equation and the Cauchy 
integral of the momentum balance equation. The Cauchy 
integral exists provided that the velocity field has the form 



v a (r,t) = —d a 9(r,t) 
m 



(54) 



where 9(r, t) is a velocity field potential. 

Starting from QHD equations we can derive an equation 
for evolution of the model function defined in terms of hy- 
drodynamic variables. Thus a macroscopic single-particle 
wave function may be defined as 



n(r, i)exp( -m9(r,t) 



(55) 
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If we differentiate this function with respect to time and 
apply QHD equations then given 

ihdMr, t) = ( -|^V 2 + T(r, t)+g \ $(r, t) \ 2 
\ 2m 



APPENDIX C 

Here we consider QHD equations for the spinning par- 
ticles and it's difference from the one's for electrically po- 
larized. 

In Euler equation contribution of magnetization M ap- 
pear instead of polarization P. Method of QHD also allow 
to obtain equation of evolution of magnetization to be 

<9 t M Q (r,t) + V^J A f (r,t) 



X V J dv"G M {v\ r")P 7 (r", t) j $(r, t), (56) 



here we use designation T(r, t) described by formula 



r(M) 



dt 



n(r',t) ' 



(57) 



The quantity T(r, t) can be considered as chemical poten- 
tial /z(r, t) and we used this fact below. The equation ob- 
tained d56]l has the form of NLSE. NLSE with DDI were 
obtained on Ref. 12311 . where the sign "-" was lost before 
integral in formulas (A6) and (A8) 

NLSE d56l ) describes collective characteristics in a sys- 
tem of many particles. This follows from the derivation 
of NLSE and from the definition of a many-particle wave 
function $(r, t) (|55]). Function T(r, t) (|57| is the contri- 
bution of the kinetic pressure and does not contain any in- 
teraction. 

If we introduce d a (r, t) as d a (r, t) = P a (r, t) /n(r, t) 
and if also we suppose that d a (r, t) is a constant d a from 
eq. ( |56l l we obtain 



ihd t $(r,t) 



2m 



V 2 + fi(r,t)+g | $(r,t) 



dPd 1 



J dr'G^{ 



r,r')|$(r',i)| 2 U( r) i) ) ( 58 ) 



where we used n(r, t) = |<3?(r, t)\ 2 . 

At temperature T equal to zero in external electric field 
E all dipoles directed parallel to external field. If E = 
E e z we have d = d e z . In this case in equation d58l) only 
one component of Green function G^ 1 of DDI is: G zz = 
(1 - 3z 2 /r 2 )/r 3 + 4vr/3 • <5(r). As z 2 /r 2 = cos 2 9 and 
in the absence of 5 function from ( |58l l we obtain the GP 
equation (TJ presented in introduction. 

The GP equation was obtained for the case of system of 
parallel dipoles. We suppose it can be approximately used 
for analyze of dipoles whose direction slowly changed in 
space. 



f&il M ' 3 (r,i 



(59) 



where J^f arises. Vanishing by thermal motion we have 
J«f = M a vP. 

For obtaining of QHD equations we started from many- 
particle Schrodinger equation 



thd t i/j s (R,t) 



((E(^-7.«?W) 



E 



u i:j - j^g: 




(R,t). (60) 



where we include short range and spin-spin interac- 
tions, and action of external magnetic field on spin. In 
Schrodinger equation d60l) we used following designations: 
7i is the gyromagnetic ratio, p ; = —iK\Ji is the opera- 
tor of momentum, fTy present short range interaction, the 
Green function of spin-spin interaction has form Gff = 

47r<5 a/} $(r 12 ) +V?V?(l/r ia ). 

For spin matrixes sf the commutation relations are 



Thereby we consider Bose particles we present here the 
evident form of the spin matrixes for particles with spin 
equal to 1 : 










For derivation of QHD equation from Schrodinger equa- 
tion we start as usual from definition of concentration of 
particles around point r of physical space: 
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where dR = Y\a=i ^ r i- ^ n this case we obtain an equations 
analogous to ([3) and © where instead polarization P(r, t) 
appear magnetization M(r, t). 
The magnetization arise in the form 



x^+(R,t) s {s^{R,t)) s . (61) 

Differentiating magnetization (I6TT ) with respect to time and 
using Schrodinger equation d60b we came to equation (l59l) . 

More details for obtaining of QHD equations for spin- 
ning particles are presented in Ref.s [23], ll69ll . 

We can perform derivation of NLSE for spinning par- 
ticles, in the result we have obtained an equation analo- 
gous to (l56l l. In the case parallel dipoles we have following 
equation, where the difference between Green function of 
spin-spin EDM interactions is accounted 

ihdMr, t) = I -t^V 2 + h(t, t)+g\ $(r, t) | 2 



/ dv' 1 J^/ m^-^f'»r.l)^ j<T>(r.M 



r — r 



8^ 

y 



(62) 

Thereby, feature of spin-spin interaction give us additional 
term in GP equation for spinning particles. This term 
caused by long-range interaction, but it has form analogous 
to SRI. 
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